Decoherence and entanglement in tripartite systems: Cavity QED 
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We propose an experimental setup, feasible with present day technology, involving two high- 
quality-factor cavities, one Ramsey zone and a two-level atom which interacts with them. The 
dynamics in the cavities is modeled by a dissipative dispersive Jaynes-Cummings model. Analytical 
expresions for the evolution superoperator in each part of the proposed setup and for the concurrence 
between the atom and the first cavity are given, and the pairwise entanglement of the subsystems 
is calculated numerically. We consider initial separable states given by the tensor product of a 
superposition of the two atomic relevant states and coherent states for the fields. We show that for 
realistic model parameters (including dissipation), when the dissipation rates are smaller than the 
dispersive frequencies (the ratio between the squared vacuum Rabi frequencies and the detunings), 
the final state of the cavity fields remains entangled. 



PACS numbers: 03.67.Bg,42.50.Pq,32.80.-t 



I. INTRODUCTION 

Entanglement [1-3], a characteristic feature of quan- 
tum mechanics, plays a key role in many quantum infor- 
mation tasks [4]. In order to implement those tasks, it 
is important to devise schemes for the creation of entan- 
gled states and to observe their dynamical characteristics 
in the presence of enviromental effects. Cavity quantum 
electrodynamics (QED) systems, the cleanest open quan- 
tum systems known so far, are interesting from the point 
of view of transfering quantum information from matter 
to radiation, from flying to static systems, and viceversa. 
Moreover, cavity QED setups [5, 6] have been used in 
many fundamental experiments in quantum mechanics, 
such as the observation of decoherence [7] , the generation 
of atomic EPR pairs [8] and the observation of quantum 
jumps [9]. 

This paper studies the entanglement dynamics under 
dissipation in a tripartite system, simple enough to allow 
for a quasi-analytical solution but rich enough to display 
a variety of dynamical characteristics, which can be im- 
plemented experimentally with present day technology. 
Figure 1 shows the system, which is composed by sev- 
eral cavities crossed by a two-level atom. Similar propos- 
als with multiple cavities, attempting to approach the 
Heisenbcrg limit for the estimation of an atomic phase 
shift [10] or producing field-field entangled states [11], 
can be found in the literature. In our system, the first 
and the last cavities are high quality cavities and the 
second one is a Ramsey zone. We calculate the degree of 
entanglement of the three pairs of subsystems, which are 
shown to behave as qubits for appropriate initial condi- 
tions, and discuss the role of the initial photon number 
and dissipation rate in the entanglement properties as 



measured by concurrence. 

This paper has four sections: the first, section II, es- 
tablishes the mathematical model that describes the pro- 
posed setup; the second (section III) presents the evolu- 
tion superoperator for the model; in the third section 
the state of the system when the atom traverses the first 
cavity and the atom-field concurrence are given; finally, 
section V describes the behavior of pairwise concurrence 
obtained numerically by using typical experimental val- 



II. MATHEMATICAL MODEL 

We consider an apparatus composed of three cavities 
crossed by a two-level atom, as depicted in fig. 1. As 
it is clear from the picture, the atom undergoes several 
dynamical processes: in the first step the atom interacts 
with a high quality dissipative cavity containing a coher- 
ent state, when the atom leaves this first cavity it will 
accumulate a phase due to free motion until the Ramsey 
zone is reached, where its state is rotated and, immedi- 
ately thereafter a second free fiight, the atom interacts 
with the final high-Q dissipative cavity prepared in a dif- 
ferent coherent state. Although the cavities contain, in 
principle, infinite modes of the electromagnetic field, we 
consider only one mode per cavity field, almost resonant 
with the relevant atomic transition. It is important to 
realize that we are assuming that the center-of-mass co- 
ordinate of the atom is treated as if it were a classical 
variable, which is possible when the velocity of the atom 
is large enough not to be affected by the changes of the 
potential energy caused by the interaction with the fields. 

We assume that the whole setup is cooled to very low 
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FIG. 1. (Color online) Scheme of the proposed experimental 
setup. 



temperatures (of the order of IK in this case) which al- 
lows us to model cavity field decay as interaction with a 
zero-temperature bath. We suppose that the atom-field 
detuning, in each cavity, is large enough to be in the dis- 
persive regime. We can describe the electromagnetic field 
of the Ramsey zone as a classical one, in its interaction 
with the atom [12]. Taking into account these consid- 
erations we can write the equation of motion of p, the 
atom-fields state, as 



at in 



(1) 



where Cd, which models the dissipative effects, is written 



field in the interval [t2,i3]- The coupling constant be- 
tween the atom and the i-th cavity field is Mli, where the 
vacuum Rabi frequency ^li is proportional to the prod- 
uct of the atomic dipole with the electric field having the 
root mean square amplitude of the vacuum in the cavity. 
The probability of detecting n photons in the i-th cav- 
ity, Pm — {a\ai\p\a\ai) ^ is usually employed to specify 
the conditions under which the dispersive approximation 
is valid. We define Ni, the relevant number of excita- 
tions in the i-th cavity, as the maximum number n, for 
which the probability P„ is still important. If the de- 
tuning of the i-th cavity field, — Ua ^ Cji, satisfies 
the inequality |Ai| Vli^/Ni + 1, we are in a region of 
parameters where the dispersive approximation is valid, 
and the atom-field interaction can be described by the 
hamiltonian [13] 



H, 



huj, ( (a|a, + 1) |e) (e| - a'^a, \g) {g\ 



(5) 



with uji — ilf/Ai. The Ramsey zone, which rotates the 
internal state of the atom, is modeled by the Hamiltonian 



Hr — e 2 
where the Hamiltonian 



'HRoe-^''^^'Wt„tdt)^ (6) 



7i (2aip{t)al - p{t)a\ai - a\aip{t)^ 



(2) 



The rate at which photons in the i-th cavity are lost, 
27i, is proportional to the cavity frequency, uii, and in- 
versely proportional to the quality factor of the cavity, 
7i cx Cbi/Qi. Since we assume that the dissipative pro- 
cesses of the atom are much slower than those of the 
field, we ignore these processes in the description of the 
dynamics of the whole system. The Hamiltonian H com- 
prises two parts, which we designate by i?o and Hi. The 
first part of the Hamiltonian, the free Hamiltonian iJo, 
is given by 



Ho = —a. 



t 1 

a,; a,; + - 



(3) 



where Wa is the transition frequency between the rele- 
vant states of the atom, \g) (ground) and |e) (excited), 
and CTz is the operator \e) {e\ — \g) {g\. The creation a| 
and annihilation operators of the cavity fields, satisfy 



the usual bosonic commutation relations 



ai,a 



= St 



— u — [a^ , ajj. The second part of the Hamil- 
tonian is the interaction Hamiltonian 



and Itti, a," 
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Hj = HnWt„,tAt) + HRWt,,fAt) + Hi2Wt,,t,{t), (4) 

where the function Wt^_t^(t) is equal to 1 if i belongs to 
the interval [ta,ib], and zero otherwise. Hence, the atom 
interacts with the first (second) cavity field in the in- 
terval of time 1^0: ^i] ([^4:^5]) and with the Ramsey-zone 



Hro = ( ^^CTz -I- hQiil^oWx ] , 



depends on flji the vacuum Rabi frequency of the Ram- 
sey zone and |^o|i the squared root of the mean number 
of photons of the coherent state sustained in the Ramsey 
zone. The operator ax is given by |e) {g\ + \g) (e|. 



III. DYNAMICS 

In order to solve the dynamics of the whole system in 
the interval [to, t^] it is convenient to find solutions in the 
intervals [ti, ti+i], for i = 0, 1, 2, 3, 4. For example, in the 
first interval we can write 

p{t)^^o{t,to)p{to)Wt„,tdt), 

where the evolution superoperator Uo{t, to) acts on p{to), 
the initial state of the total system. The function W 
stresses that the last equation is valid only in the interval 
[to, ti]. The state of the system at the end of one interval 
is the initial state for the next interval. Thus, we can 
write, in general 



pit)^i^kit,tk)pitk)Wt„t,^,it), 



0,1,. ..,4. 



Employing Lie algebraic techniques [14], as shown in 
appendix, we cast the evolution superoperator in the 
[tk,tk+i] interval as 



Uk{t,tk) = Uk{t,tk) ■ Ul{t,tk) Uokit), 



(7) 
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where Uk{t,tk) is the unitary operator which solves the 
Hamiltonian dynamics 

ihjUk{t,tk) = {Ho + Hi{t)) Uk{t,tk), Uk{tk,tk) = I, 

in that time interval, and Uok (t) is a superoperator given 
below. Wo employ the dot-convention which establishes 
how (super)operators act on other (super)operators: the 
(super)operator to the right of the superoperator with 
dots will occupy the place of the dots. For example, 

A ■ B{p) = ApB, A-{B-C)=AB-C. 

The dissipative part of the evolution superoperator reads 



= 1,2 



where the operator J-ik{i) (i = 1, 2; fc = 0, 1, 2, 3, 4) 

^ik{t) = — 7, —■ 7 N (8) 

27^ + iLOik [az <Tz ) 

depends on the parameters wife = LUiSkfi-, W2ft = W2(5fe,4. 
The field operators Mi,Vi and Ji are given explicitly by 



Mi = alai-, Vi^-ajai, and = • oj. (9) 

The unitary contribution to the evolution superoperator 
can be cast as 

Uk{t,tk) = e-'^'-'''^^"°'y'^Uk{t,tk) 

where Uk{t-tk) is the hamiltonian evolution operator in 
the interaction picture. Hence, t/i = / = C/3 is the iden- 
tity in the intervals of free fiight; 

Uk{t, tk) = e-'(*-*'=)[^"/^='-l/^, fc = 0, 4, 

while the atom crosses the high-Q cavities and 

U2{t,t2) = e-'^^*-*^)^^ exp [-iftRl^olit - t2)a4 , 

while the atom crosses the Ramsey zone, assumed to be 
resonant with the atomic transition. 



IV. ENTANGLEMENT AND DECOHERENCE 
ON THE FIRST CAVITY 

The initial state of the two-cavities-atom system that 
we consider, which is relatively easy to be prepared with 
current experimental techniques, is a separable one: the 
atom in a superposition of its two relevant states and the 
fields in coherent states, p{to) = |V'(*o)) {'4'{to)\, where 

|V'(io)) = ^(|e) + e-^/^|<?))®|a)®|/3), 



and I a) is a coherent state for the first (second) 

cavity. While the atom crosses the first cavity, the state 
of the system is 

pi{t) = ^{\e,al){e,al\ + xit)\e,al){g,al\ (10) 
+x* {t) \g, a?) (e, a? | + \g, a?) {g, af |} ® (^1 1 , 



where 



a. 



ae 



a? = ae 



-iCj\t—iu)\t—^\t 
-icjit-\-icJit—yit 



(11) 



_ aa-iC02t-j2t 



and 



x{t) = exp 



71 



71 + icoi 



-i-^ + laff^t) - + LUa)t , (12) 

_(1 _ e-2(7l-Hiu;i)f) _ _ e-27lt)_ 



We have chosen to = 0. 

Since the atom-first cavity field density operator be- 
longs to a (g) subspace of the total available Hilbert 
space, it corresponds to a couple of qubits, whose entan- 
glement can be quantified by means of Wootters concur- 
rence [2]. In our case 



C{M = Kt)|<^l-|(af|af)|2. (13) 

Clearly, one of the factors on which concurrence depends, 

is a function of the nonorthogonality of the states ja^) 
and I Off), the more different the states the larger the con- 
currence. Although the other presents a more complex 
dependence, its main contribution is of the form 



-|ap(l- 



-271* 



that is, it principally depends on the initial number of 

photons and the product of the decoherence rate and the 
time spent in the cavity. In fig. 2 we show the time 
evolution of the concurrence between the atom and the 
first cavity for different values of initial mean number of 
photons and different dissipation rates, the time interval 
was chosen between to = and ti = 1000 iis. As the 
initial average photon number is increased, the maximum 
entanglement lasts for longer times in a lossless cavity. 
However, at the same time, for a lossy cavity, at larger 
values of a, atom and cavity field separate faster. 

A graphic way to see this is by looking at the phase 
space representation of the field [17]. The physics of this 
part of the process seems easy to understand: without 
dissipation, the key element for entanglement, the inter- 
action, is proportional to the field intensity. However, for 
coherent states, the role of the interaction is to split the 
field into a component correlated to the cxc;itcd atomic 
state whose label rotates at angular velocity uii + uii and 
a component correlated to the ground atomic state whose 
label rotates at angular velocity u)i — oji (See eq. (11)). 




Phase space evolution 

mean value of a, 



FIG. 2. (Color online) Atom-first cavity field concurrence 
while the atom crosses the first cavity, as a function of the 
initial mean number of photons, (a) a =0.5, (b) a =1, (c) 
a =2 and for different values of the dissipation rate. In the 
second column a blow-up of the first two-oscillations is shown. 
For to = and ti = 1000 /^s. 



It is easier to concentrate on the angular velocity differ- 
ences, and to plot the labels in a phase space picture (fig. 
3). It is clear that every 27r/ (2a;i) seconds both labels co- 
incide, hence the state becomes separable. This process 
remains unaffected in the presence of decoherence. From 
eq. (13) it is clear that entanglement decreases as the 
overlap increases, or, in other words, entanglement de- 
creases with the phase space distance between the field 
states correlated with the excited and ground states, dis- 
tance which is given by the length of the vertical line 
that join those states. For larger mean number of pho- 
tons, the phase space distance is large most of the time, 
except near the separation times; for small mean number 
of photons this distance never grows appreciable, neither 
does the entanglement. 

In the presence of dissipation the characteristic time 
for decoherence is inversely proportional to the average 
photon number. We have already pointed out that a 
factor which decays with the characteristic time of deco- 
herence enters in the expression for concurrence. In fig. 3 
one observes that the part of the field entangled with |e) 
(|g)) rotates clockwise (counterclockwise) and its radius 
shrinks with dissipation. The assymetry of the entan- 
glement curves shown in the fig. 2 can be understood 
in terms of the phase space distance, as was discussed 
above, and marked on fig. 3 with vertical lines. 

As the atom crosses the second cavity, the expressions 
for the density operator calculated analytically become 
much more cumbersome and are omitted. Therefore, con- 
currence is calculated only numerically. From here on we 
make a more numerical approach. Taking into account 
typical values of current experiments [7] we consider the 
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FIG. 3. (Color online) Phase-space evolution of oi\ and oi\ in 
a rotating frame with frequency iDi . The initial mean number 
of photons is Q = 1. Dashed line shows the dissipationless 
case and solid lines represents the dissipation one with 71 = 
1.25 kHz and dispersive frequency uj\ = 6.25 kHz. Vertical 
distance of the trajectories at same time is proportional to 
entanglement. The line AA'{BB') corresponds to a time of 
larger (smaller) entanglement. 



following fixed parameters: the angular dispersive fre- 
quencies uji = UJ2 = 6.25 kHz^ the atomic transition 
frequency and the Ramsey zone frequency Wa = ^^i?, = 
51.1 X 10*^ kHz, the detuning |Ai| = IA2I = 100 kHz and 
Rabi frequencies Q,i — Q.2 — 25 kHz. The flight times 
through the cavities are also fixed: the atom spents 30 /is 
in each cavity and 10 us between cavities. The dissipa- 
tion rates for the cavities 71 and 72 are varied according 
to 71 = gwi and 72 = qui2, where g,q = 0, 0.05, 0.5, 1. 
In order not to break the dispersive approximation, we 
can not increase too much the mean number of photons. 
To estimate the validity of the dispersive approximation 
we use |Ai| ^ Sly^jap -I- 1. For the values of the param- 
eters chosen above, we are constrained to coherent states 
satisfying |ap '-^ n ^ 15. We thus choose, for the mean 
number of photons, the values a, /3 = 0.5, 1, 2. 



V. ENTANGLEMENT AND DECOHERENCE 
ON THE SECOND CAVITY 

We begin by considering the first cavity to be ideal and 
the second to range between ideal (q=0) and strongly dis- 
sipative (q=l). In fig. 4 we plot the concurrences for all 
three pairs of subsystems. For ideal cavities, and initial 
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FIG. 4. (Color online) Atom-fieldl (a) atom-field2 (b) fieldl- 
field2 (c) concurrence for ideal first cavity and different dissi- 
pation rates for the second one. 



coherent states with a — 0.5 — (3, the final state contains 
quantum correlations between all pair of subsystems. If 
the second cavity is strongly dissipative only the atom 
and the first field show relevant quantum correlations. 

If we consider the opposite case, that of a ideal second 
cavity, and first cavity ranging between ideal and strongly 
dissipative, entanglement behavior is complementary to 
the previous case (fig. 5): only the atom-second cavity 
field are important for a strongly dissipative first cavity. 
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FIG. 5. (Color online) Atom-fieldl (a) atom-field2 (b) fieldl- 
field2 (c) concurrence for ideal second cavity and different 
dissipation rates for the first one. 

We conclude that the most interesting physical effects 
from the point of view of multipartite entanglement hap- 
pen for decoherence rates of about one order of magni- 
tude smaller than the dispersive frequency. This conclu- 
sion holds even if we change the mean number of photons 
in the cavities, as shown in fig. 6. All these results shows 
that the local independent interactions between the atom 
and the cavity fields has as consequence the generation 
of non-local field- field entanglement. 

As noticed before, increasing the mean photon num- 
ber can lead to smaller atom-cavity (quantum) correla- 
tions (see fig. 7). In fact, atom-first field concurrence 
decreases as the initial mean number of photons in sec- 
ond cavity increases from one to four. Also, if we begin 
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FIG. 6. (Color online) Fieldf-Field2 concurrence. Maximum 
values are obtained for (nearly) ideal cavities and minimum 
values for strongly dissipative cavities. 



with coherent states with mean photon numbers larger 
than one, even for not so large values of the decoherence 
rate, we observe sudden death of the atom-first field en- 
tanglement (see fig. 7, first row). As shown in figures 4-5 
and 7, the entanglement between all pairs of subsystems 
can not increase simultaneously, suggesting some kind of 
entanglement conservation law, related to the Coffman- 
Kundu-Wootters (CKW) monogamy inequality [18]. 




FIG. 7. (Color online) Atom-fieldl, atom-field2 and fieldl- 
field2 concurrence for an ideal first cavity, a = 0.5 and differ- 
ent values of q and /3. 

It is interesting to notice that as the initial mean num- 
ber of photons increase so does the field-field entangle- 
ment, for (nearly) ideal cavities, and short times. As 
time goes on, this rule does not hold anymore, and there 
seems to be an optimal value of the mean number of pho- 
tons to produce the higher field-field entanglement for a 
given interaction time (see fig. 7, for g=0 and q=0, 0.5). 

In summary, we have considered a realistic cavity QED 
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setup, assumed to be operated on the dispersive regime, 
and have solved analyticaUy its dynamics, in the sense of 
finding the evolution superoperator. We have shown that 
for typical experimental values of the cavity dissipation 
rates, smaller than the dispersive frequencies, and times 
of the order of 100 fis, not only the final state generally 
exhibits atom- field and field- field entanglement, but that 
the effect of the dissipation is not trivial, in the sense that 
the pairwise concurrences are not monotonic functions of 
the dissipation rates. 
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Appendix: Evolution superoperator 

It is easy to show that the evolution superoperator 
Uk{t,tk) in the interval [tk,tk+i], k = 1,2,3, given by 

Texp/^.(/:, + l[i^,.] + f [i^.M,.]), 

where T is the time-ordering operator, can be written as 
(7) due to the following commutation relations 

[Ho;Cd]=0=[-Ho,CdI 
[HR{t);CD] = 0=[-HR{t),jCD]. 

In the intervals [to,*!] and [t4,t5] the evolution superop- 
erator 

Uk{t,tk) = e'^^'-"'\ fc = 0,4, 
can be written as 

U{t,tk) = )e^'^ \ (A.l) 



with 

= [Ha + Hf + '^ uiik [alci 

i 



O'z ■ + ■ O'z 

2 

(^z) 



+ \e) {e\,-], 



(Mi + Vi) , 



where ojik = ijJiSk,o, ijJ2k = W2(5fc,4 and the operators 
M.i-,'Pi,Ji are defined in (9). Taking into account that 
the atomic operator commute with the field operators, 

we write 



^^ocik{Mi+Vi)+l3ikJi 



where the commutation relations 

[Mi + Vi, Jj] = -2JjSij, i,j = 1, 2, 

and disentangling operators techniques [14] were used. It 
is clear that the coefficients Uik contain a term propor- 
tional to 7i and one proportional to oJik, which can be 
separated. We can write therefore 



,aik(n=0)(Mi+Vi) 



(A.2) 



gaiS,(a;ifc=0)(Xi+7'i)g|^^(exp(2a.fc)-l)Ji 



Grouping the first term of (A.l) with the terms above 
for zero dissipation (7j = 0) we recover the hamiltonian 
dynamics. On the other hand the last line of (A.2) cor- 
responds to the purely dissipative dynamics (with zero 
Hamiltonian). After employing the explicit expressions 
for aik,Pik, the evolution superoperator can be written 
as (7). 
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